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<^> ' Abstract 

oo ; 

Q\ . This paper analyses the Noether symmetries and the corresponding 

' conservation laws for Chern-Simons Lagrangians in dimension d = 3. 

, In particular, we find an expression for the superpotential of Chern- 

O 1 Simons gravity. As a by-product the general discussion of superpoten- 

1} , tials for 3rd order natural and quasi-natural theories is also given. 
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^ . 1 Introduction 

The interplay between symmetries, conservation laws and variational prin- 
ciples for Lagrangian dynamical systems have been under consideration for 
a long time, starting from the celebrated work by E. Noether. Our aim in 
this paper is to calculate conserved quantities for Chern-Simons gravity in 
dimension d = 3. Chern-Simons theories exhibit in fact many interesting 
and important properties. First of all, they are based on secondary charac- 
teristic classes discovered in ||] and there is some hope that they will give 
new topological invariants for knots and three-manifolds. Moreover, they 
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appeared in physics as natural mass terms for gauge theories and for gravity 
in dimension three, and after quantization they lead to a quantized coupling 
constant as well as a mass [Q]. They have also found applications to a large 
variety of physical problems Q, among which we quote anyons and quan- 
tum Hall effect; when the gauge group is the inhomogeneous Lorentz group 
then such a theory is equivalent to standard gravity. Chern-Simons gauge 
theory is also an example of a topological field theory M. It should also be 
remarked that, as it has been recently shown, Chern-Simons Lagrangians 
are the only obstructions to an equivariant inverse problem of the calculus 
of variations in d = 3 |5|. Finally, the Chern-Simons term is related to the 
anomaly cancellation problem in d = 2 conformal field theories ||. 

Our paper is organized as follows. For our later convenience, we shall 
start from a short, elementary and self-contained exposition of a variational 
approach to the problem of conservation laws in field theories. A more 
systematic approach can be found e.g. in @, |, @, [fij>0» while for a rigorous 
mathematical presentation [] we refer the reader to [l^, 13, 14, [l^]. In section 



3 we collect some material from our earlier papers concerning superpotentials 
and in particular we shall list the explicit formulae which are necessary to 
calculate them in all theories of order at most three. Section 4 deals with 
Chern-Simons gauge Lagrangians and we discuss how general formulae for 
superpotentials apply in this case. Some results of section 4 are finally used 
in section 5 to calculate the energy-momentum complex for Chern-Simons 
gravity. 



2 Preliminaries 

Let us consider a field theory over a spacetime manifold M of arbitrary 
dimension d, describing the behaviour of an arbitrary field by means of a 
Lagrangian L = L((p, d x (j), . . .) of arbitrary order. The action functional will 
be 

I(<t>) = I Ldx 
Jd 

where the integration is performed over a domain D of the manifold M of 
independent variables x = (x M ); the field configuration (i.e., the set of all 
dependent variables) is represented here by (f). Variation (i.e. a functional 
derivative) 51 of / with respect to an arbitrary (infinitesimal) variation 8(f) 

1 It involves a variational calculus on jet bundles. 

2 For simplicity we drop an internal field index, e.g. (f> A 
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of tp leads to the variation 5L of the Lagrangian itself. This, in turn, can be 
rearranged into two parts according to the "first variation formula" (] 

5L= 6 ^ + d^ (1) 

The first term consists of the Euler-Lagrange expression (field equations): 

5L dL dL dL 

YI = ^I- + ^Ktt— ) - • • • 



and therefore vanishes on shell, i.e. when the field tp satisfies the equations of 
motion. The second part (boundary term) is a divergence of g^ = g^(tp, Sep), 
where 

depf, dtp dtp^ 

This boundary term is usally neglected since it does not contribute to the 
equations of motion (and it is killed on the boundary dD after imposing 
suitable boundary conditions on 5t/)). However, this second contribution is 
physically important: in fact, it does contribute to conservation laws. 

Recall that a variation, say 5*tp, is called an (infinitesimal) symmetry of 
the action / if the corresponding variation <5*L can be written as a divergence 
(i.e., 5*L = d^T^) without using the equations of motion. Of course, this 
depends heavily on the transformation properties of L: for example, if the 
Lagrangian transforms as a scalar (i.e. L{tp) = L(tf) + 5*tp)), then one has 
5*L = 0. 

Therefore, for the variation implemented by a symmetry transformation 
(5*0, equation (1) can be rewritten under the following form 

5 ±5^ = -d,{ e »-T») (3) 

where g* = g^^tp^^tp). A Noether current then arises 

E» = E»(<l>,8.<f>) = <g-TK (4) 
which is conserved on shell. One writes 

d^E* = mod — 



Here the usual Einstein summation convention over repeated indices is adopted. 
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or, equivalently, 

dpE» « 

and calls it a weak conservation law. It follows then from the Poincare 
lemma that (at least locally) there exists on shell a skew-symmetric quantity 
U^ v = —U v ^, called a superpotential, such that 

i.e. E^ differs from the divergence d^U^ u by a quantity which vanishes on 
shell. Notice that the superpotential so defined depends explicitly on field 
equations. However, if a formula on U^ u can be "analytically" prolonged 
to all field configurations^ then a quantity e M = d u U^ u defined for all field 
configurations is automatically conserved also "off shell", i.e. 

V = 

independently of the dynamics. In this case we call this a strong conservation 
law. 

3 Second Theorem of Noether and Superpoten- 
tials 

In the present paper we shall be interested in so-called local symmetries 
(and second Noether 's theorem). This corresponds to the situation when 
the variation 5*(j) depends on some (finite) number of arbitrary functions 
9 = {9 l (x)} i = 1, . . . , N over spacetime M together with their derivatives 
up to some (finite) order, say s. In other words one has 

8*4, = 5 e = + tfdpO 1 + ... + ^- ps d pi . . . d Ps e l (5) 

The order s of the highest derivative involved in this field transformation 
law is called the geometric order of the transformation. If the Lagrangian 
is of order k, i.e. it involves up to k derivatives of the fields, the integer 
r = s + k is then called the total order of the theory |l2|] . 

Supposing now |^ that the total order of the theory is at most 3, then 
both sides of equation (3) become differential operators acting on 9, where 

e» = e»(6) = tfe 1 + tfdpQ 1 + t^dpdj 1 (6) 

4 This is always the case for so-called local symmetries (see below). 
5 We slightly extend here the formalism developed in where proofs of our statements 
can also be found. 
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^S e <j> = W(9) = WiP + wedpO* + wTdpdrP (7) 

Notice that in (6) and (7), in order to preserve the uniqueness of the de- 
composition, one should assume the following symmetry conditions for the 
highest coefficients |: 

Replacing now (6) and (7) into equation (3) and calculating the divergence 
d^E^ explicitly, one gets an identity which holds true for arbitrary func- 
tions B % . Therefore the total coefficients resulting in front of 6 l , d a 9 l , etc., 
should vanish independently. This produces, in fact, a differential inden- 
tity involving the Euler-Lagrange expression, so-called generalized Bianchi 
identity [13], which in this case reads as follows 



Wl - d v w\ + d v d pW \ p = (8) 

The coefficients of W as given by (7) vanish on shell because of the Euler- 
Lagrange equations, so that dpE^ ~ holds. It appears, in this case, that 
the Noether current admits a representation of the form 

E f " = E f " + dpU pp ^iP + eP (9) 



where E 1 * is the "reduced" Noether current, which vanishes on shell (sec fl3|| ). 
The reduced current E is unique, while the superpotential is not, since (9) 
remains unchanged if U^ u is redefined by the addition of the divergence of 
any skew-symmetric quantity Y^ vp = yt'w], it can be found that equation 
(9) has an explicit solution given by 

& = « - dpw^y + wfdpd' (io) 
uw = (t M + a^iwy + tfdvP (ii) 

where tf pv = (4/3)t[ wl ". This gives explicit formulae for a superpotential in 
the case of theories of total order at most 3. It should be noticed that the 
"lowest" coefficients tf do not enter this explicit expression (11); therefore, 
this suggests that superpotentials are, in a suitable sense, "algebraically 
simpler" then the corresponding Noether currents. 

We should also remark that formula (11) seems to provide a "canonical" 
expression for the superpotential; but, in general, this is not true. In fact, 

6 We shall use the notation A 1 ^ = ^(A^ - A vti ) for the skew-symmetrisation and 
A (^) = U A iw + for the symmetrisation. 
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equation (11) explicitly depends on the local coordinate system we have 
chosen as well as the transformation properties of the coefficients appearing 
in (6). Nevertheless, as we shall see below, such a canonical expression will 
be possible under some additional assumption. 

As a general class of examples including theories of gravity (which in 
fact we shall need in the sequel) let us assume that the field is a field 
of geometric objects over M |l2| , i.e. we assume that under coordinate 
changes in M the field <fi transforms properly through laws which depend 
only on a finite number of the partial derivatives of the coordinate change 
itself Q . We assume also that the Lagrangian is natural (or, equivalently, 
reparametrization invariant) , in the sense that the lift of all diffeomorphisms 
of M transform L as a scalar density of weight 1. This means that, at the 
infinitesimal level, one has 

8*(j) = £(</> (12) 

and 

5*L = £{.L = d a {i a L) (13) 

where £ = £, a d a is an arbitrary vectorfield on M (i.e., an infinitesimal dif- 
feomorphism) and denotes the Lie derivative along £. 
As a consequence, equation (3) can be written as follows: 

dpE* = 8^(1), Stf) - eL] = - 6 ^C^ 

where = £7**(£) is the appropriate Noether current evaluated along the 
infinitesimal symmetry given by £, called now energy-momentum flow. It is 
known that the Lie derivative of any geometric object can be expressed by 
means of covariant derivatives with respect to an arbitrary symmetric linear 
connection on M, In particular, we may expand the Lie derivative C^4> as a 
linear combination of £ and its covariant derivatives, i.e. we write 

= * a e° + $£v P r + • • • (14) 

Suppose again ^ that the total order of the theory is at most 3. This implies 
global and covariant decompositions 

E H _ = + T ^ Vp e Q + T^VpVaC (15) 

7 Tensors or linear connections are geometric objects, while arbitrary gauge fields are 
not. 

8 This includes all standard metric theories of gravity (k — 2, s = 1) as well as so- 
called first-order formalism (or Palatini), with k = 1, s — 2 (Yang-Mills theories have 
k = 1, s = 1). 
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s ^8^ = w(0 = w a e + wpy P e + wp°v p v a e (ie) 



According to the general techniques developed in [13] we get thence covariant 
versions of equations (9), (10) and (11), i.e. 

W„ - V V W£ + V v V„Wp = (17) 

as a generalized Bianchi identity; 

E" = (W£ - V p W^)e + W£ p V p e (18) 

as a reduced energy flow; and finally 

jjw = ( T M + v v j%\rt)g* + f^ v V v i a (19) 

where fg pv = (4/2)Tt p]u (recall that T^ [H = 0). Equation (19) gives 
the required formula for a canonical superpotential in the case of natural 
theories of total order at most 3 (see [[HJ). 

As an example one can consider the superpotential obtained from the 
Einstein-Hilbert (purely metric) gravitational Lagrangian (here notation is 
standard) 

L EH {g,dg,d 2 g) = R„ v 

which is given by 

U£ u H (0 = \detg\ 1 H^e-^e) 

and it is known as the Komar superpotential 0. It has been recently shown 
that the Komar expression is "universal" in the following sense: for a large 
class of non- linear gravitational Lagrangians (in the first-order "Palatini" 
formalism), the superpotential does not depend on the Lagrangian and is in 
fact equal to the Komar superpotential |l6|, |l7]]. 

In some mathematical literature currents admitting superpotentials are 
called trivial since they lead to strong conservation laws, which hold true 
irrespectively of the form of the equations of motion [18]. Such currents are 
also somehow "trivial" from a co-cohomological viewpoint (in the sense the 
they define trivial co-homology classes in de Rahm cohomology, being on 

9 We refere the reader to e.g. [| for more exhausted discussion of the Komar 
superpotential. 

10 It has been also found that in these cases the universality property holds for the 
Einstein equations as well. 
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shell exact forms). Nevertheless, we remark that they produce physically 
relevant quantities, like e.g. charges, masses and so on. 

The advantage of using superpotentials is twofold. First, as we have 
mentioned before, they are in general (algebraically) simpler expressions 
then the corresponding Noether currents. Moreover, in order to calculate 
the flux of any conserved quantity they allow, via Stokes' theorem, to reduce 
calculations to the following: 

JT, JdT, 

Here E is an hypersurface in M of codimension 1 and <9£ denotes its bound- 
ary (which has codimension 2). In the case of d = 4 spacetime this simply 
means that the flux of a conserved quantity through a 3-dimensional portion 
of spacetime can be calculated as the integral of the superpotential along 
the 2-dimensional boundary S£, which is a surface. 



4 Chern— Simons Gauge Theory 

To start with, we discuss another physically interesting example, i.e. the 
Chern-Simons theory considered as a gauge theory over M = M 3 (see also 
[§, |l0|). The Chern-Simons three- form written in terms of a connection 
and its curvature reads as follows 

L^ = tr(flAw--a?AwAw) (20) 

where u> = A^dx^ is a connection one-form, is a matrix- valued R gauge 
potential, Q = du> + u> A lj = ^F^dx 11 A dx u is its curvature two-form, 
F,ui; = d^A u — dyA^ + [A M , A u ] is the gauge field strength and tr denotes 
the matrix trace. 

The corresponding Lagrangian density Lq S = Lcsdx 1 /\dx 2 Adx 3 written 
in terms of more physically relevant quantities is 

Lcs = \e^tx{¥^A p - ^A M A,A p ) (21) 

where is the totally skew-symmetric Levi-Civita symbol in d = 3. 

11 In fact, A M belongs to the Lie algebra of some gauge group G, which for simplicity 
we assume here to be represented in terms of (complex) matrices. 
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In this case the first variation formula (1) reads as 

5L CS = e^tr(F pu 5A p ) - d p [s^tr( AJA p )} (22) 

and this gives F pu = as the equation of motion, which means $7 = 0, i.e. 
that only flat connections are solutions of the equations of motion. 

As a symmetry one can consider any infinitesimal gauge transformation 

5 x A p = D pX (23) 

where x is an arbitrary (matrix- valued) function and D p = 9L + [A^,-] 
denotes the covariant derivative with respect to the connection u. It is 
known from transformation properties of the Lagrangian Lcs ( see e -g- Si) 
that the following holds: 

5 X L CS = d p [e^tr(A u d pX )] (24) 

Replacing now (22), (23) and (24) into equation (3) we get 

e^tr{Y^D p x) = d^E cs { X ) (25) 

where 

E CS { X ) = 2e^tr(A,A pX ) + 2e^tr(A v d pX ) = (26) 

-2^tr(A,A pX ) + 2e^tr(A u D pX ) 

is the expression for the Noether current stemming from the gauge invariance 
of lcs- Therefore, indentifying the corresponding coefficients in (26) with 
those in (7), we obtain with the help of (11) the following superpotential 
(see also M): 

Uj&ix) = 2e^tT{A vX ) (27) 

The left hand side of (25), after expanding it as a differential operator acting 
on x, gives rise (see also equation (8)) to the generalized Bianchi identity 

e^D p F pu = (28) 

which, in dimension d = 3, is known to be equivalent to the standard Bianchi 
identity for the gauge field F pu , i.e. 

DpF pi/ + D p F vp + D u F pp = 
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A similar analysis can be performed for the diffeomorphism invariance 
of Lcs\ in this case one has 0: 

5 ( A p = £ € A„ = eOaA^ + A a d^ a (29) 

and 

5(.Lc S = Ct-Lcs = d a (t a L cs ) (30) 

Now we have: 

E»(0 = eLcs + e^tr{A u d a A p )e + e^tr{A v A a )d p e (31) 
from which it follows via (11) 

U^s(0 = e^tr{A v A a )e (32) 

since only the last term in (31) contributes into (32). 

The generalized Bianchi identity takes thence the form 

e^tr(D p F^A a - F^F ap ) = (33) 

in which the first term vanishes owing to (28), so that (33) reduces in fact 
to the funny algebraic identity 

£ ^tr(F^F ap ) = (34) 

which holds true for any gauge field in dimension three! 

For a comparison, we recall the reader that the gauge and diffeomorphism 
invariance of Yang-Mills Lagrangians in four dimensions lead instead to 

U^ M (x) = tr(F^ X ) (35) 

UZAO = tr(F^A a )r (36) 
Some autors [20, 21] prefer to deal with so called improved diffeomor- 
phisms 

<%A P = e a F ap = <%A P - D^Aa) (37) 
which differ from the Lie derivative (29) only by a gauge transformation 
X = £ a A a . Direct calculations show that 

U£ P S (0 = U£ p s (0 - U£UX) = £ ^tr(A,A a )r (38) 
i.e., it differs from (32) by a sign. Similarly, one finds from (35-37) 

U£ P M (0 = (39) 

It means that there are no, in the Yang-Mills theory, non-trivial charges 
corresponding to these transformations. 

12 Under diffeomorphisms the gauge potential A M behaves as a one-form and Lcs as a 
scalar density of weight one. 
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5 Chern— Simons Gravity 



We now proceed to discuss a less studied case of Chern-Simons gravity con- 
sidered as a purely metric theory over M 3 . Therefore, the only dynamical 
variable is a Riemannian (or pseudo-Riemannian) metric g = g^ v . The 
Lagrangian density is 

Lcsg = ie^"tr(R^r p - ^r M i\,r p ) 

1 2 

_ -e fjjUp (Rp fll/ T ap - -r^r^rap) (40) 

where, using the notation introduced in the previous section, we have now 
u = T^dx^, is the matrix-valued one-form of Riemannian connection, 
= ^Tl^dxt 1 A dx v is the Riemann curvature two-form and Yi^ u = d^Y u — 
d u T ^ + [r^, T u ] is the Riemann curvature tensor written in matrix notation. 
For the sake of completeness we also recall two basic formulae of Riemannian 
geometry 

= \g aa (dpg, a + d^p - d a g Pll ) (41) 

Rppv = — dv^pn + r^r^ — r^r^ (42) 

as well as the following indentity 

Rp^v = d^Rpv - $vRpn + 9(3uR c f l - gp^RZ + ^Ridp^Z ~ 9pvd%) ( 43 ) 

which holds true only in dimension d = 3 (see e.g. |l9[| ) f^. Here, R^ u = R^ au 
denotes the Ricci tensor and R = i?" is the Ricci scalar of g. 

Notice that Lcsg depends on the metric g only through the Riemannian 
connection (41). The relation (22) remains still valid 

5L C SG = £^tr(R^r p ) - d^ v Ptx{T v 8T p )] 

but it does not correspond now to the first variation formula (1), since T p 
is no longer a dynamical variable in our theory. Accordingly, one has to 

13 Here T M = Fp^ and R M „ = Ji«„„ are represented as 3 x 3 matrices and the upper 
index is the row index. 

14 This means, in fact, that in three dimensions the two notions of Ricci flatness and 
flatness coincide: i.e., that gravity is trivial in dimension d = 3. This is why a coupling 
with Chern-Simons term enriches the theory. 

15 We shall be frequently using g for raising and lowering indices. 
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replace 5T p in the first term by means of the well known " Palatini formula" 

ST$p = \f° ( Vptgpa + V p 5g a p - V a 8g Pp ) (44) 

Making use of equation (43) together with the contracted Bianchi identity 
V a R* = \V p R one finds 

e^Ri pv bV% = -2e^V^5g ap + d p (2e^R«5g ap ) (45) 
With the help of the tricky identity which holds true for any quantity S pu 

2s^S% = ^(S^-S^) (46) 

the symmetric part of the first term of the right hand side in (45) can be 
converted into the following form 

-2e^V p (R: - ^R)Sg ap 
Therefore, the first variation formula reads now as 

5Lcsg = -2C a P5g ap + d p \e^(2R a v bg ap - tr(T u 5T p ))} (47) 

where 

C ap = e^ p V p (R« - ^8£R) (48) 

is the so-called Cotton tensor ^] , the vanishing of which gives, in the case 
of Lcsg, the Euler-Lagrange equations of motion. It is also known (see e.g. 
[§, [H]]) that the Cotton tensor is symmetric, traceless (C° = 0), divergence- 
free (V a Cp = 0) and it vanishes if and only if the metric g is conformally 
flat (= conformal to a flat one). 

As a symmetry transformation, consider then a 1-parameter group of 
difieomorphisms generated by the vectorfield £ = £ a d a . In this case one can 
use the well known expressions 

<%g = £t9ap = Va£p + Vpia (49) 

and 

SjTp = C^ ip = CRLp + VaV p f (50) 

16 In fact, C ap is a tensor denity of weight 1. 
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In a Riemannian case we are dealing with, the last formula is, of course, 
equivalent to the Palatini one, if one replaces 5g in (44) by (49). 

Thus, our task is now to calculate S^Lqsg directly from its transforma- 
tion properties. Notice that the transformation rule for the linear connection 

dx x „ n dx a dx 13 dx x d 2 x p 



N _ 



^ <9a^ a/3 dx^ dx' 4 da^ dx v dx p 

,dx x p dx a dx x d dx p dx 13 
~^dx~P a/3 dx^ + d^dx^^dx^^'dx^ 1 
is the composition of a vector- and of a pure gauge transformation. In other 
words one has 

f^jr^ + j^j- 1 )— (si) 

where J = is the Jacobian matrix. Therefore, at the infinitesimal level 
(x a = x a — £ a ), S^Lcsg is the sum of two terms of the type (24) and 
(13). More precisely we have obtained 

S^Lcsg = d P [FLcsg + e^tr(T u d pX )] (52) 

where x = g§a is an infinitesimal version of the inverse Jacobian matrix. 
We see that the Lagrangian Lcsg is not natural, but it is in some sense 
quasi-natural, i.e. it differs from a natural one by a total derivative. 
In a similar way one deduces from (51) that the following holds 

S^T p = V p x + ed a T p + T a d p C (53) 

as the sum of (23) and (29). Here V p x denotes the "formal" covariant 
differential of %i formally treated as a tensor of rank (1, 1), i.e. we set 
^ pX = 9 p x + [r p ,x]. Of course, it can be directly checked that equations 
(50) and (53) are equivalent. 

Inserting now all necessary expressions in equation (3) we obtain the 
variation formula under the form 

2C a P^g ap = -dpE£ SG (0 (54) 

where 

E csg(0 = PLcsg + e^f {tr (T u d a Tp)C - 2R T Mg pa + 5°g Ta )V (r e+ 
17 In fact, we have to substitute \ l— * ffa i n equation (24). 
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+tr(r„[r p , x]) + tr(r,r a )5 P r + 2tr(r„0„x)} (55) 

is the energy-momentum complex for the Chern-Simons gravitational La- 
grangian (40) . We are now in position to obtain the corresponding superpo- 
tential U G P SG . Since superpotentials are additive quantities one can calculate 
the contributions from all terms of (55) separately. Moreover, the first two 
terms are of the lowest order and do not contribute to the superpotential. 

To obtain the contribution from the third term of (55) one will use the 
covariant formula (19). First we rewrite this term under the form T^ a 'V a ^ a , 
where T pa = —2e pup R^,(5^g pa + d^g^a). Then, by making use of the iden- 
tity (46), tI m<t ' can be converted into the explicit skew-symmetric form 
e pau (3R ua - Rg ua ). This gives rise to: 

3U CSG = ^ P (R9u a ~ 3R„ a )C (56) 
Similarly, contributions coming from the fourth and the fifth terms are 
easy to calculate directly by using (11) and (46): 

4 U£ P SG + 5 U^ SG = e^{[T a ,T u r a + tr(T u T a )}e (57) 

where [IV, = T^T^ - T^ u T^ a , which due to the symmetry of T^ a 
reduces (57) to the form: 

4 u£ p SG + 5 u^ SG = e^{r^rt)t a (58) 

The sixth and last term in (55) is of the second order. We write it 
as t ppa d p d a £, a , with t%> a = 2e^ p Tal- One needs the shorthand = 
4/3ta^ a (see (11)); a direct calculation with the help of another tricky 
identity 

gives then rise to = 2e^r^+(2/3)e WCT rf ^ Since % U£ P SG = 
t^d^, one gets: 

6 U^a = l^ Up (d,T% - Sfyrt)^ + \er">{W w - 8^)8^ (60) 

We can now formulate our main result: the superpotential U G P SG (£) ap- 
pears as the sum: 

U£ P sg(0 = 3 U^ SG + 4 U^ SG + 5 f/^ G + *U PP SG (61) 

with addenda given by (56), (58) and (60). 

The left hand side of (54) taking into account (49) has the form 4C^V Q £*. 
This, of course, leads via (17) to a divergence-free property of the Cotton 
tensor as a generalized Bianchi identity. 
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